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Abstract 


We study the distribution function of the 6-quark for inclusive semilep- 
tonic B decays in the framework of the parton model. We compare the func¬ 
tional behavior of the fragmentation inspired function and the Gaussian AC¬ 
CMM function, and discuss their validity. Using the HQET parameters, we 
can critically test these models. 
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I. INTRODUCTION 


Studies of inclusive semileptonic decays of B mesons have drawn our interest as a testing 
ground for the motion of 6-quarks inside hadrons. We have the differential decay width for 
the process B ^ Xqlu 
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where is the leptonic tensor and the hadronic tensor is dehned as 

= \Y. j - q - p,){B\J>‘^\X)(X\r\B). (2) 

spins \ '^) ^ 

Since the B meson is heavy, the momentum transfer of the decay is larger than the typical 
energy scale of hadronic binding, ~ ^qcd, in rnost regions of phase space. This enables 
us to formulate the inclusive semileptonic decay process in an analogous ma nn er to the 
well-known method of the deep inelastic scattering (DIS). Light-cone dynamics dominates 
the decay process and the hadronic form factors show the remarkable feature of scaling in 
the absence of the relevant corrections. So the parton picture should work well, where the 
modeling of the decay process is possible by summing the partonic decay processes weighted 
by the distribution of the 6-quark momentum inside the B meson. Bareiss and Paschos 
noticed this point and they formulated the inclusive decays along the lines of analyzing DIS 
in the parton picture. The general expression for the hadronic tensor in terms of invariants 
is given by 
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Following the parton picture of Feynman |Q] , we let the probability of hnding a 6-quark in a 
B meson carrying a fraction x of the meson momentum in the inhnite momentum frame be 
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the distribution function f{x). The decay width is then obtained by incoherently summing 
the partonic decay widths weighted by f{x) in that frame. Thus we immediately hnd that 

g ■ Pfl) = I if{x+) + fix-)) 

W^ 2 (g^ q-PB) = ix+fix+) - x-fix-)) 

W"3(g^ q-PB) = ifix+) - fix-)) 

W^4,5,6(g^g■pJ = o (4) 


after neglecting the lepton masses where 

go ± y'q^ + m] 


x± = 




( 5 ) 


where niq is the hnal state quark mass in 6 — qeu. This approach has been reformulated on 
the light-cone in Ref. [^]-|^ and holds in any frame. The QCD radiative correction of the 
parton model was studied in Ref. P . The idea of the parton model was also applied to the 
exclusive semileptonic decays of B meson in Ref. |^. 

Recently the heavy quark effective theory (HQET) has been applied to inclusive semilep¬ 
tonic decays of B mesons and brought great progress im- With the expansion for the 
decay rate in inverse power of 6-quark mass, a systematic study is possible and one obtains 
more information about the hadronic matrix elements than before. Perhaps one of the most 
important results of the HQET analysis for the inclusive decays is the absence of the con¬ 
tribution of the l/ruft correction terms. With the help of the l/rub expansion, we can relate 
the moments of the charged lepton energy spectrum to the matrix elements for the operator 
expansion of the current product. 

The l/mfo structure of the inclusive parton model has recently been analyzed in Ref. []^ . 
When the parton model approach motivated by DIS is applied to the analysis of B decays, 
the modeling boils down to introducing the distribution function of Eq. (4), which describe 
the 6-quark’s motion inside the B meson. In other words, in the study of the inclusive 
semileptonic B decays, the problem of which model we choose corresponds to the problem 
of choice of the distribution function. In general two models are known which predict the 
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lepton energy spectrum of inclusive decays of B meson. Altarelli et al. (ACCMM) 
treat the fe-quark and the spectator quark as quasi-free particles with a Gaussian spectrum 
of Fermi momentum in the B rest frame. Bareiss and Paschos 00 notice that the b- 
quark distribution inside a B meson in the inhnite momentum frame (IMF) is related to 
the fragmentation of a 6-quark into a B meson by crossing symmetry of the spectator, 
and so their functional behaviors are similar to each other. Thus they chose Peterson’s 


fragmentation function as the distribution function of 6-quark, which is the default 
choice in Lund Monte Carlo programs. Both models get results at the hadronic level after 
convolution of the quark level spectrum. The difference between the two models mainly 
comes from the weight function which is corresponding to the distribution function f{x) of 
6-quark momentum in the parton model approach. 

In this paper, we study both models from the viewpoint of the parton model picture 
to compare their distribution functions with each other. For explicit comparison of the 
distribution functions of both models, we derive the distribution functions described by the 
same definition in the same reference frame. Moreover, with the help of the l/m^ study 
suggested in Ref. we attempt to test the models by extracting the HQET parameters, 
the effective mass of light degrees of freedom A and the average kinetic energy of 6-quark 
We show that known bounds on A and /i^ constrain the model parameters and may 
even rule out the model itself. In section 2, we review the parton model approach hrst. In 
this approach Peterson’s function and its improved form by Lee and Kim 0 is taken as the 
distribution function. In the ACCMM model, the distribution of the 6-quark is determined 
by the Gaussian distribution for the Fermi motion of the spectator quark in the B meson 
rest frame. In section 3, we derive the distribution function of the ACCMM model in the 
inhnite momentum frame from the Gaussian function, which allows a consistent comparison 
of the two models. The behavior of the distribution functions is investigated and the heavy 
quark effective theory (HQET) parameters are obtained in terms of the parameters of each 
model. Our discussions and conclusions are summarized in section 4. 


4 













II. THE PARTON MODEL 


We first review the parton model approach applied to analyze inclusive semileptonic B 
meson decays. As described in the previous section, the distribution function f{x) is defined 
as the probability of hnding a 6-quark in a i? meson carrying a fraction x of the meson 
momentum in the inhnite momentum frame. Thus the decay probability is determined by 
the momentum of the 6-quark and has the value of f{x)dx. And we write the Lorentz 
invariant decay width as follows 

Eb dT(B —> XqCf) = J dx f{x) Et dT(b qeu), (6) 


with the relation = xpg. Since B mesons are heavy, light-cone dynamics dictates the 
inclusive semileptonic decay analogous to that of DIS. The universal distribution function 
of QCD is interpreted as the structure function of the 6-quark in the parton model, which 
determines the distribution of the light-cone projection of the 6-quark momentum inside 
the B meson. We directly derive the triple differential decay rate as 

dr _Gl\Vqk\^ qo-Ei 


dEidq^dqo + 


[x+f{x+){2Ei - mgxE) + {x+ ^ a;_)] 


(7) 


rat 


where x± is dehned in the previous section. Jin and Paschos argued that the contribution 
of f{x-) is expected to be relatively small in the kinematic region we are interested in, and 
can be neglected. This is a consequence of the fact that the light-cone distribution function 
is sharply peaked around x ~ mi,/nig. 

For a more quantitative analysis, we have to introduce a method analogous to the HQET. 
In this parton model approach, information on the nonperturbative dynamics inside the B 
meson is encoded in the distribution function. Therefore, it is crucial to determine precisely 
the form of the distribution function in this approach. Unfortunately we don’t presently 
know how to calculate the true form of it. However, we have some knowledge about the 
line shape function of QCD from the help of HQET and the operate product expansion 
(OPE). Since the mass of the 6-quark does not appear in this parton model framework, we 


5 





have to define it within this framework in order to express the model in terms of the inverse 
power expansion of the 6-quark mass. With properly defined rrib, the parton model is well 
represented by an 1/mb expansion |^. We define 

mb = {x)mg, ( 8 ) 


where (x) = J dx xf{x). With this definition as the starting point, we can relate the parton 
model to HQET and compute the effective mass of the light degrees of freedom in a i? 
meson, A. 

In the HQET framework, the triple differential decay rate for B Xqlu shows scaling 
behavior 0-0: 

dEidq^dqo 1927r3 A ® mb — qo ’ 

where 


ml + q^ - 2mbqo 
2A (mb - qo) 


( 10 ) 


Apart from kinematic factors, the triple differential decay rate is completely determined by 
one scaling variable Xg. The scaling behavior of QCD will be violated by perturbative and 
nonperturbative (higher twist) corrections. If we take the final state quark to be w-quark 
and so m^ = 0 hereafter, the scaling form is derived from Eq. (6) 

dr ^ GllKbP 2 24{qo-Ei){2mbEi-q^) 

dEidq^dqo 1927r3 m^ mb — qo 

where m^^ is as defined in Eq. (7). From Eqs. (8) and (10), the line shape function E{Xg) 
of QCD is related to the distribution function f{x+) of the parton model after appropriate 
variable changes 


= jFim- 


( 12 ) 


Next we intend to study their functional behavior. Following the notation of Bigi et ah 
HJim , a general form for the singular expansion of the distribution function is written as 
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where moments of the equation: 


a„ = y dxg xlF{Xg), n = 0,1, 


(14) 


Using the HQET and the OPE, we parametrize the moments with the forward scattering 
matrix elements of B meson. Especially we know that oq = 1 from 6-number conservation, 
and that Ui = 0 up to 0{l/ml), which indicates the lack of 0{l/mh) corrections. The value 
of the hrst moment Oi = 0 is related to the value of the 6-quark mass or the parameter 

A. 

The second moment of the spectrum is related to the average kinetic energy of 6-quark 
inside the B meson and usually represented by the parameter /i^, dehned by 


0-2 — 


^ ^ -{Blbn^blB) = 


(15) 


' 3A2’ 

where tTq = iDa — rubVa with the covariant derivative and Va the four-velocity of the 
B meson. After the variable changes of Eq. (11), we can calculate the moments for the 
distribution function in the framework of the parton model. Under the dehnition of mb in 
Eq. (7), we have Oi = 0. The second moment is derived as 

«2 = dx^ xlF{x^) = y dx {x- {{x'^) - {x)^). ( 16 ) 

We hrst used Peterson’s fragmentation function as the distribution function of 6-quark 
following the original work of Bareiss and Paschos [^, and obtained the value for 02 ~ 0.76, 
with the model parameter e = 0.004 and the mass of B meson = 5.3 GeV, which agrees 


well with the QCD sum rule result of Ref. |T^. As pointed out in the Ref. ||I^, however, 
this function gives rather larger A and /i^, and we need to hnd a better functional form 
or some modihcation for Peterson’s function. Jin and Paschos 0 also pointed out that the 
distribution function should be more sharply peaked than Peterson’s function. By improving 
Peterson’s arguments, Lee and Kim [Q suggested a modihed form of Peterson’s function 
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(17) 


/r(^) = 


N, 
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Using this new form, we calculate A and /i^ by varying the input parameters and e^. 

Fig. 1 shows the behavior of the distribution function with varying and e^. The 
larger the value of a^, the sharper the functional behavior. And the smaller the value of 
Cq, the sharper the function and the closer the location of peak to 1. We also show the 
parameter space of for allowed A and /i^ in Fig. 2. For the values of /i^, there is 

some theoretical controversy which we will discuss later in section 4. We take the range of 


A and /i^ obtained by the HQET and QCD sum rules in Refs. [IT,IT 


A = 0.4 ~ 0.6, 


/i; = 0.6±0.1 [16], and /i^ = 0.10 ± 0.05 


(18) 


Ball and Braun [|T^ calculated using the QCD sum rule approach and obtained = 
0.60±0.10 GeV^ for R-meson, while Neubert 0 obtained —Ai = 0.10±0.05 GeV^. We hnd 


that there exists a region in the ( 0 ^, 6 ^) plane that gives values of A and in agreement 
with those of Ball and Braun [^. If Neubert’s prediction [^ for the value of /i^ is correct, 
there is no common region in which the model parameters give the values of A and in 
Eq. (17), and we may conclude that this functional form is not appropriate for inclusive 
semileptonic B decays. 


III. THE ACCMM MODEL 


The AGGMM model considers the B meson as consisting of a 6-quark and a spectator 
quark. The spectator quark is treated as a quasi-free particle having dehnite mass m^p, and 
momentum p, while the 6-quark is treated as a virtual particle of invariant mass mb 


ml(p) = m| + mjj, - + ml^, 


(19) 


where is the B meson mass and p = |p| in the B meson rest frame. The 6-quark energy 
is given by 

8 






































Eb{p) = \/p^ + ml{p) 


(20) 


= m^- ^p^ + mip . 

The distribution of the spectator momentum is usually assumed to be a Gaussian with an 
adjustable parameter 
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$(p) = 




/ p- 


with the normalization 


poo 

/ dp = 1. 

Jo 


( 21 ) 


( 22 ) 


After convoluting the distribution of the spectator momentum with the subprocess of the 
virtual 6-quark decay in a frame where its momentum is —p, we can obtain the lepton 
energy spectrum at the hadronic level. This model has been extensively used in the analysis 
of the charged lepton energy spectrum in semileptonic decays and numerically reproduce 
the inclusive spectra well. In the light of QCD, several authors analyzed the model 
and suggested a way to reproduce the pattern of the QCD description of the heavy quark 
motion. In this paper we derive the distribution function f[x) in the IMF from the given 
Gaussian distribution in the B rest frame, ‘h(p), in order to study the ACCMM model by 
comparing with the parton model approach. 

Let’s choose a coordinate system where the B meson is at the origin in B rest frame. 
We boost the B meson along the positive z direction with the velocity /?. As is well known, 
for the B meson at rest we have the boosted energy of B meson as 


= = (23) 

where 7 = 1 /and that of 6-quark as 

E;{p) = jiE,ip)-Pp,), (24) 

where p is the magnitude of the spatial momentum of the 6-quark and pz is the 2 ; component 
of p in the 5-rest frame. Now we calculate the energy ratio between 6-quark and B meson 
in this moving frame 
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(25) 


E*, 7{E,-Pp,) 

E*b iras 

In the limit of /3 —1, 7 —> cx), we define the ratio as 


El 


{Eb - Pz) 


m. 


= X, 


(26) 


which represents the ratio of 6 -quark energy to B meson energy in the IMF. 

Because no specihc direction for p is preferred in B meson decay, the value of Pz is 
distributed from —p to +p with equal probability for a hxed value of p. Then the value of 
X can be from Xminip) to Xmaxip) with equal weight for a hxed p, where 


Xminip) = - {Eb{p) -p), 

rrig 

{p) = —{Eb{p)+p), 


Xr 


(27) 


with the weight m^/2p, because Xmaxip) — Xminip) = ‘^p/'^b- Each value of p in the B 
meson rest frame corresponds to a distribution of x in the IMF such that the probability of 
hnding a 6 -quark in a i? meson carrying a fraction x of the meson momentum in the IMF is 
obtained by intergrating over p with Gaussian weight, ^(p), 

f'Pmax ^ ‘YYl 


rp 

fix) = / 

Jo 


dp p^ ^{p) e{x- Xminip)) diXmaxip) “ x) 

Zp 

Pmax _ 777 , 

dp p2 *.(p). ^ 


(28) 


^ Pmin 


where 


Pmin 


|m|(l - xf - 7 n^p| 


(29) 


2 m^(l — x) 

is obtained by the conditions x^inip) < x < Xmaxip)- We take the value of Pmax to be 
inhnity as far as the integration is concerned, because the area of phase space for p > Pmax 
is negligible. After the integration over p, we obtain 


fix) = 


m. 


V^Pi 


■ exp 


PPf 




(1 -x) 


(30) 


4 \mg{l - x) Pp 

where p = rnfgp/p],- Using the Eq. (12) and the relation, (1 — x+) = ;^(1 — x^j), we hnd 
Eq. (30) becomes Eq. (30) of Ref. 0, the light-cone distribution function for the ACCMM 
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model derived by Bigi et al. We show the functional behavior of f{x) in Fig. 3 by varying 
the parameters and Note that the shape of the function f{x) depends on two 

parameters rrisp and p^. The peak position is determined by msp- 


^peak 


rrig - rrisp 


(31) 


while the Fermi motion parameter p^ is related to the width of the distribution /(ac), i.e. 
the larger the value of the broader the distribution f{x). One can easily check the 
normalization of f{x) 


f{x)dx = 


rPma 


dp p^<h(p) = 1. 


(32) 


The hrst moment of f{x) is directly related to the average 6-quark energy in the B rest 
frame, 


(x) ^ f' dx xf(x) = r™' dp pMp)^ ^ 

Jo Jo m„ m„ 


(33) 


because rrif, = {x)m^ in the parton model approach we get the relation = {Eb)B-rest, 
which reminds us of Csaki and Randall’s dehnition of the 6-quark mas^ for the ACCMM 
model with respect to the HQET. 

We thus have derived a distribution function f{x) of Eq. (28) for the parton model 
approach from a Gaussian distribution in the B rest frame. Following the method discussed 
in section 2, we can calculate the moments of the line shape function and the parameters A 
and for given model parameters. First we calculate A : 


A = mg-mb = m^- {Eh)B-rest (34) 

roc - - 

= rrig — J dp p^<F(p)(mg — m^p) 

= i^p"^ + m‘^^p)B-rest- 


^ Csaki and Randall showed that the predictions of the ACCMM model agree well with those 
of the heavy quark effective theory with a definition of the 6-quark mass, mf, = {Ei,{p))B-rest- 
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From Eq. (14), is given as 


And using 


we get 






ix^) = [ dx x^f{x) 
Jo 


9 ii^b)B—rest q (P )B—rest)i 




Ftt ~ ^{{Eb) B-rest i^b) B-rest) A {p ) B-rest- 


( 35 ) 


(36) 


(37) 


Finally with the following relations, 

{El) B-rest - {Eb)l_rest = {P^)B-rest + rn% - {^ld^^^^p)l-rest (38) 

and 

nCO Q 

{p"^) B-rest = dp p^^{p) = -pI, (39) 

we obtain 


pi = Spl + 3m% - 3(y'p2 + 


( 40 ) 


We show the parameter space of {Pp,msp) for given A and pi in Fig. 4. Values of (p^, m^p) 


in the shaded region A satisfy the condition of Eq. (17) given in the Ref. |^, while values 


in the region B satisfy the condition of Eq. (17) given in the Ref. |^. We can determine 
the model parameters Pp and rusp, once we know the correct value of the HQET parameter 




IV. DISCUSSION AND CONCLUSIONS 

We have studied two phenomenological models (QCD parton model and ACCMM model) 
for the analysis of inclusive semileptonic B decays in the context of the parton model ap¬ 
proach. In this framework, the core of modeling is the choice of a distribution function of 
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the 6 -quark inside a B meson. The light-cone dominance enables us to describe the triple 
differential decay rate in terms of a single function, the light-cone distribution function up 
to kinematic factors, as the form given in Eqs. (9) and (11). The distribution function in 
ACCMM model comes from the Gaussian distribution of the 6 -quark momentum at the B 
rest frame. Bigi et al. [Q obtained the light-cone distribution function for the ACCMM 
model by reading off it from the photon energy spectrum of i? —>■ Xg'y decays. We calculate 
it here by counting the probability of finding a 6 -quark with the boosted energy ratio x 
between 6 -quark and B meson into the infinite momentum frame within the framework of 
the parton model. We find that our result of Eq. (30) exactly reproduce the light-cone 
distribution function given by Bigi et ai, Eq. (30) of Ref. [§]. This agreement is originated 
from the fact that the ACCMM model treats the 6 -quark as a quasi-free particle at the de¬ 
cay rate level by classical treatment as we do not consider the perturbative QCD corrections 
in this paper. 

Even though we cannot calculate the distribution function from QCD completely, some 
model independent information can be obtained with the help the HQET. When we expand 
this function with infinite number of moments, we cannot know all the moments while the 
HQET give some informations for a few moments. We note that since this is an singular 
expansion, we cannot construct a well-behaved function with only a few moments. In spite of 
that, the information is very useful for giving some constraints on phenomenological models. 
Investigating the distribution function, two moments of the function, corresponding to the 
mean value and the width, can be related to two QCD parameters, A (or mb) and /i^. The 
mean value of the function strongly depends on A and very weakly on the /i^, while the 
width on the /i^. Actually the fact that ai = 0 define the quark mass, equivalent to A and 
the next moment 02 is directly related to the /r^. With these constraints, we obtain the 
possible regions of model parameters for each model, as shown in Figs. 2 and 4. 

Recently, it has been an important subject to obtain an accurate value of the kinetic 
energy, (= — Ai), of the heavy quark inside B meson in connection with the heavy quark 


effective theory (HQET). Neubert derived |T^ the field-theory version of the virial theorem 
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within the HQET framework and obtained the result ~ 0.10 ±0.05 GeV^, which is much 
smaller than the QCD sum rule calculations of Ball and Braun and rather comparable to 
the earlier QCD sum rule result using a less sophisticated approach in Ref. ||2^. However, it 
should be noted that Refs. and differ in the choice of the 3-point correlation functions 
used to estimate the matrix elements of interest. The difference in the numerical values 
obtained in these two calculations is understood in terms of the contributions of excited 
states, which in principle must be subtracted in any QCD sum rule analysis. In practice, this 
subtraction can only be done approximately. Therefore, the numerical differences between 
Refs. and indicate the limited accuracy of the QCD sum rule approach. Bigi et al. 
m derived an inequality between the expectation value of the kinetic energy of the heavy 
quark inside the hadron and that of the chromomagnetic operator, (p^) > |(My2-Mp2), 
which gives /i^ > 0.36 GeV^ for B meson system. However, Kapustin et al. m showed 
later that this lower bound could be significantly weakened by higher order perturbative 
corrections. By means of a QCD relativistic potential model, the value of /i^ was found to 
be quite large as 0.44 ~ 0.46 GeV^ in Ref. |2^. Besides the above theoretical calculations of 
Gremm et al. extracted the average kinetic energy by comparing the prediction of 
the HQET with the shape of the inclusive B —*■ Xlv lepton energy spectrum for Ei > 1.5 
GeV, in order to avoid the contamination from the secondary leptons of cascade decays of 
6 —^ c —>■ slu. They obtained —Ai = 0.19 ± 0.10 GeV^. Combining the experimental data 
on the inclusive decays of D Xeu, B ^ Xev and B Xtu, Ligeti et al. a derived 
the bound of /i^ < 0.63 GeV^ if A > 0.240 GeV, or < 0.10 GeV^ if A > 0.500 MeV. Li et 
al. [0 obtained the value of —Ai centered at 0.71 GeV^ from the analysis of the inclusive 


radiative decay B —> within the perturbative QCD framework. Related with the 

comparison of various theoretical calculations of we note that Ref. emphasizes that 
one has to be careful when comparing the values of —Ai obtained using different theoretical 
methods. The mixing of the operator for the heavy quark kinetic energy with the identity 
operator occurs at the two-loop order and the parameter Ai requires a non-perturbative 
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subtraction. 


Note also that the functional behavior of Refs. 0 and m is different from each other 
while both functions satisfy the constraints. Even though two moments, ai, 02 , which 
correspond to the mean values and the widths of the distribution functions, agree with 
each other, their full behaviors do not necessarily agree. For the two models studied in 
this paper, we hnd that the modihed Peterson’s function in Fig. 1 has a rather longer tail 
than the distribution function derived from the Gaussian function of the Fermi momentum 
dehned in the ACCMM model in Fig. 3. We cannot say which model is better, at least 
at this stage. Both models can give proper values for the 6-quark mass and the width 
of the charged lepton energy spectrum, both theoretically and phenomenologically. For 
further tests, we must £t them to the experimental data more precisely after more data is 
accumulated. For this purpose the photon energy spectrum of the inclusive radiative decay 
B Xs'j would be more useful. The photon energy spectrum of this process is known 
to be more sensitive to the behavior of the distribution function than the charged lepton 
energy spectrum of the inclusive semileptonic decays [^,^, since the charged lepton energy 
spectrum is contaminated by kinematics. 
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FIGURES 

FIG. 1. (a) Functional behavior of the improved Peterson’s function with varying = 0.05, 

0.085, 0.1 for fixed = 0.004; (b) the same function with varying = 0.002, 0.004, 0.006 for 
fixed Oq = 0.085. 

FIG. 2. The model parameter space in the {a^, e^) plane for the QGD sum rule results 
for A and Values of (a^jCg) in the shaded region satisfy the condition 0.4 < A < 0.6 and 
0.5 <nl< 0.7. 

FIG. 3. (a) Functional behavior of the distribution function in the infinite momentum frame 

derived from the Gaussian distribution of spectator momenta with varying = 0.3, 0.5, 0.7 GeV 
for fixed nisp = 0.1 GeV; (b) the same function with varying rusp = 0.0, 0.2, 0.4 GeV for fixed Pp 
= 0.5 GeV 

FIG. 4. The model parameter space in the (pp, rUsp) plane for the QGD sum rule results for 
A and Pp. Values of {pp,msp) in the shaded region A satisfy the condition 0.4 < A < 0.6 and 
0.5 < Pp < 0.7. Values of {pp,msp) in the shaded region B satisfy the condition 0.4 < A < 0.6 and 
0.05 <pI< 0.15. 
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